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Because of the controversial X(3872) meson's very close proximity to 

the D U D threshold, this charmonium-like resonance is often considered 
a meson-meson molecule. However, a molecular wave function must be 

essentially of a meson-meson type, viz. D U D in this case, with no other 
significant components. We address this issue by employing a simple two- 
channel Schrodinger model, in which the J PC = 1 ++ cc and D°D* chan- 
nels can communicate via the 3 -Po mechanism, mimicked by string breaking 
at a sharp distance a. Thus, wave functions and their probabilities are 
computed, for different bound-state pole positions approaching the D°D 
threshold from below. We conclude that at the PDG X(3872) mass and 
for reasonable values of a, viz. 2.0-3.0 GeV -1 , the cc component remains 
quite substantial and certainly not negligible, despite accounting for only 
about 6-10% of the total wave-function probability, owing to the naturally 

long tail of the D D component. 

PACS numbers: 14.40.Pq, 14.40.Lb, 12.40.Yx, 11.80.Gw 

The X(3872) charmonium-like meson is by now a very well established reso- 
nance [1]. It was first observed in 2003, by the Belle Collaboration [2], in the 
decay — > K^it + it~ J/ip, with significance in excess of 10<7. Since then, 
it has been confirmed by several collaborations, viz. Belle, BaBar, CDF, 
DO, CLEO, and, more recently, by LHCb (see the 2012 PDG pQ listings 
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for references). The PDG summary table lists the X(3872) as an isoscalar 
state with positive C-parity, from the observed 7 J/ip decay, but unknown J 
and P, having an average mass m = 3871.68 ± 0.17 MeV/c 2 and a width 
r < 1.2 MeV/c 2 . The two most likely J PC assignments are 1 ++ and 2 '", 

while the observed hadronic decay modes are p°J/ip, ojJ/ip, D°D* , and 

D°D°7r . Henceforth, we shall denote D°D*° simply by D°D*°. 

Meson spectroscopists have been puzzled by the X(3872), because of 
its low mass as compared to predictions of conventional quark models, as 
well as its remarkable proximity to the D°D*° threshold, being "bound" by 
only 0.15 MeV pp. This has led to a plethora of model descriptions of the 
X(3872), viz. as a cc state, meson-meson (MM) molecule, tetraquark, or 
hybrid meson. For a number of reviews on the many different approaches 
and the experimental situation, see [3]. Recently, we have described [4 J 
the X(3872) as a regular but "unquenched" 1 ++ ( 3 -F*i) charmonium meson, 
whose physical mass is dynamically shifted about 100 MeV downwards from 
the bare 2 3 Pi cc state due to its strong coupling to the S- wave D°D*° 
and D^D*^ channels, besides several other OZI-allowed and OZI-forbidden 
(p°J/ip, ujJ/ip) channels. Thus, the observed hadronic X(3872) properties 
were well reproduced [4]. 

Nevertheless, the closeness of the X(3872) to the D°D*° threshold seems 
to favour a molecular interpretation [5]. In the latter paper, it is stated 
that, whatever the original mechanism generating the resonance, a near- 
threshold bound state will always have a molecular structure. This implies 
that the MM component of the wave function, i.e., D°D*°, should be the 
only relevant one. Here, we shall study this issue in a simplified, coordinate- 
space version of the model employed in [4], restricted to the most important 
channels, viz. cc and D°D*°. Note that, even if the X(3872) is essentially 
a molecule, it will mix with cc states having the same quantum numbers. 

Now we turn to the two-channel model used in [6], with parameters 
adjusted for the X(3872). Consider a coupled qq-M\M<i system, with the 
qq pair confined through a harmonic-oscilator (HO) potential, whereas the 
two mesons M\,M2 are free. The correponding 2x2 radial Schrodinger 
equation is given by Eq. ([!]), with the Hamiltonians ^ and Here, p c j 
is the reduced mass in either channel, m q = mq the constituent quark mass, 
l c , If the orbital angular momenta, and u the HO frequency: 

1 / d 2 l c (lc + l)\ 1 2 2 
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Note that we use here relativistic definitions for the MM reduced mass /if 
and relative momentum k, even below threshold contrary to [6], though 
this is practically immaterial for the X(3872). At some "string-breaking" 
distance a, transitions between the two channels are described by an off- 
diagonal point-like potential with strength g 

V = -^—5(r - a) . (4) 



Continuity and twice integrating Eqs. (|lj|3|) yields the boundary conditions 
u'Jr \ a) — u'(r I a) + —Uf(a) = u' f (r f a) — u'Ar J, a)-\ — —u c (a) = 0, (5) 
u c (r f a) = u c (r J. a) and Uf(r f a) = -uj(r J, a) . (6) 



A general solution to this problem is given by Eqs. ([7]) and Q for the 
confined and the MM state, respectively. The two-component function 
u{r) = (u c (r),Uf(r)) is related to the radial wave function as u{r) = rR(r): 



u f (r) 



Uc(r) 



A f Ji(kr) 



A c F c (r) r <a 
B c G c (r) r> a 



(7) 



B 



J lf (kr)k 2l f+± cot (6 lf (E)) - N lf (kr) 



r < a 



r > a 



(8) 



Now, F c (r) vanishes at the origin and G c (r) falls off exponentially for r 
oo. Defining then z = fitor 2 and 



E-2m c lc + 3/2 



we get 



F(r) 



T(Z + 3/2) 



2uj 



z (l+l)/2 e -z/2 



-u,l + 3/2,z) 



G(r) = -] 2 T{-v)rz l / 2 e- z/2 ^{-v, I + 3/2, z) . 



(9) 



(10) 



11) 



Here, the functions <p and ij) are the confluent hypergeometric functions of 
first and second kind, respectively, and the T function acts as a normalising 
function. The functions J and N in Eq. (|8|) are defined in terms of the 
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spherical Bessel and Neumann functions j, n, i.e., Ji( kr) = k~ l rji(kr) and 
Ni(kr) = k l+1 rni(kr). Prom the boundary conditions (5 6) and the explicit 
wave-function expressions in Eqs. (|7|8|), we obtain 



G' c (a)Fc{a) - F' c {a)G c (a) = lJ lf (ka)F c (a)% , 
J' lf (ka)N lf (ka) - J lf (ka)N[ f (ka) = *%Ji f (ka)F c (a)% 
Using next the Wronskian relations 



W(F c {a),G c (a)) 
W(N lf (ka),J lf (ka)) 
yields 



Jim [F c (r)G' c (r) - F' c {r)G c {r)] = 1 , 
N lf (kr) J' lf (kr) - N[ (kr) J t/ (kr) 



(13) 



lim 

r— ¥a 



A f B 



He 



A r B r 



and 



9 - 2 ^J? f (ka)F?(a) 
a ji c 



Be 
A, 



(14) 
(15) 



Finally, with the expression for the MM scattering wave function Uf(r) 
(second line in Eq. (pi)), the final result for cot 8\ f {E) is obtained, reading 



cot (5 lf (E)) 



g 2 ^kjf(ka)F c (a)G c (a) 

He 



1 n lf (ka) 
3i f (ka) 



(16) 



Now, in the present X(3872) model, there is only one scattering channel, viz. 
for the D°D*° system. Thus, poles in the 5"-matrix, which represent possible 
resonances, bound states, or virtual states, are given by the simple relation 
cot 5i f (E) = i. On the other hand, the solutions to the two-component 
radial wave function (7j8) are then fully determined by relations (14) and 
(15), up to an overall normalisation constant. 

Next we apply this formalism to the coupled cc-D°D*° system. In the 
confined channel, the cc system is in a 2 3 Pi state, and so l c = 1, whereas 
the D°D*° channel has If = 0. In Table [j] we give the fixed parameters of 
the model, with the HO frequency lv and the constituent charm mass as in 

Table 1. Fixed model parameters [7] and D°D*° threshold. 



Parameter 


0J 


m c 


m D o m D *o 


m D + 7720*0 


Value (MeV) 


190 


1562 


1864.84 | 2006.97 


3871.81 



[7], being unaltered ever since. However, the radial quantum number v in 
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Table 2. Bound and virtual states near the D°D*° threshold. 



a (GeV- 1 ) 


g 


Pole (MeV) 


Type of bound state 


2.0 


1.133 


3871.68 


virtual 


2.0 


1.150 


3871.81 


virtual 


2.0 


1.153 


3871.81 


real 


2.0 


1.170 


3871.68 


real 


3.0 


2.097 


3871.68 


virtual 


3.0 


2.144 


3871.81 


virtual 


3.0 


2.150 


3871.81 


real 


3.0 


2.199 


3871.68 


real 



Eq. (|9| varies as a function of the energy, and therefore will generally be 
non-integer, becoming even complex for resonance poles. The parameter 
that determines such variations is the coupling g. In the uncoupled case, 
i.e., for 5 = 0, one recovers the bare 3 Pi HO spectrum, with energies (3599 + 
2nu) MeV (n = 0, 1,2, . . .). The only other free parameter is the string- 
breaking distance a. Now we try to find S-matrix poles as a function of 
the coupling g and for two reasonable values of a, viz. 2.0 and 3.0 GeV -1 
(~0.4 and 0.6 fm). Searching near the D°D*° threshold, a dynamical pole 
is found, either on the first Riemann sheet, corresponding to a bound state, 
or on the second one, which represents a virtual state (see Ref. [1], second 
paper). These results are presented in Tableland Fig. [Tj 



Fig. 1. Dynamical real (solid) and virtual (dashed) pole trajectories for a = 2.0 
GeV -1 (left) and a = 3.0 GeV" 1 (right). The arrows indicate pole movement for 
increasing g. The PDG [Tj X(3872) mass is labelled by *. Also see Table [2j 
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Fig. 2. Radial wave-functions for E = 3871.68 MeV and g = 1.170, g = 2.199 for 
a = 2.0 GeV" 1 (left) and a = 3.0 GeV" 1 (right). Also see Tabled 




r(GeV- 1 ) r(GeV-') 



Note that the dynamical pole arises from the D°D*° continuum and 
is not connected to the bare 2 3 Pi cc state at 3979 MeV, contrary to the 
situation in [3] (first paper). For our study here, this is of little consequence. 

Finally, we depict the normalised two-component wave-function R(r) in 
Fig.[2j evaluated for the PDG PQ X(3872) mass of 3872.68 MeV. One clearly 
sees the P-wave behaviour of the cc component, whereas the D°D*° is in 
an S-w&ve. Moreover, the cc admixture is certainly not negligible, despite 
the low total probablities of 6.13% and 10.20%, for a = 2 GeV^ 1 and a = 3 
GeV -1 , respectively, which are logical because of the very long tail of the 
D°D*° component; also see [8]. Soon we will publish more detailed work. 
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